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Abstract. We conjecture the true rate of growth of the maximum size of the Ricmann zeta- 
function and other L-functions. We support our conjecture using arguments from random 
matrix theory, conjectures for moments of L-functions, and also by assuming a random model 
for the primes. 



1. Introduction and statement of results 

A fundamental problem in analytic number theory is to calculate the maximum size of L- 
functions in the critical strip. For example, the importance of the Lindelof Hypothesis, which is 
a consequence of the Riemann Hypothesis, is that it provides at least a crude estimate for the 
maximum in the case of the Riemann zeta-function. In this paper we use a variety of methods 
to conjecture the true rate of growth. 

Consider first the Riemann zeta-function, which is a prototypical L-function. The Lindelof 
Hypothesis asserts that for every e > 0, C(§ + = 0(t e ) (here we assume t is positive). Under 
the Ricmann Hypothesis, one can show that 



C(i + 1 ;) = o(cx P (c-^)) (i.D 



log log t 

for some constant C (see Theorem 14.14A of ^5], for example). Several results of the form 



have also been established. Assuming the Riemann Hypothesis (RH), Montgomery ^1 showed 
C > 1/20. Balasubramanian and Ramachandra [2] improved the constant C and removed the 
assumption of RH. Soundararajan ^Sj further improved the estimate to C > 1, and he has also 
obtained similar results for the central values of a family of L-functions. In fact, Soundararajan's 
calculations show that the proportion of t for which £ (| +it) is this big is quite large, suggesting 
that it may get bigger still. Numerical calculations of Kotnik indicate that C > 2 and 
perhaps C can be much larger. 

We are interested in finding out which of equations or l|1.2f> is closer to the truth. This is 
part of a class of problems that has recently come to be known as the "1 or 2?" question, where 
one has an O-result and an ^-result which, suitably interpreted, differ by a factor of 2. In the 
case here, the unknown factor is the power of logt in the exponential. The calculations in this 
paper support the view that "1" is the correct answer in this case, and we make the following 
conjecture: 
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Conjecture A. 

max |C(| + it) | = cxp ((1 + (l))^ilo g rioglogr) . (1.3) 

Similar arguments to those presented for £(A +it) work for S(t), the error term in the number 
of zeros of the zeta-function up to height t, and lead to 

Conjecture B. 

l imsU p S{t) = _L (1.4) 

t— >oo yj log t log log £ TTy/2 

Much recent progress in understanding analytic properties of L-functions has come from 
the idea of a "family" of L-functions with an associated symmetry type. The idea is that to 
a collection of L-functions, with appropriate natural conditions, one can associate a classical 
compact group: unitary, symplectic, or orthogonal. One expects the analytic properties of the 
L-functions to be largely governed only by the symmetry type. Here we apply this philosophy 
to conjecture the maximal size of the critical values of L-functions. 

A family T of L-functions is partially ordered by the "conductor" c(F) for F € T . Our 
calculations assume that #{L € T : c(F) < D} w D. Straightforward modifications can 
handle the case in which the family grows like D A for any A > 0. 

For a more detailed discussion of families of L-functions, see (However, note that [S] 

introduces a refined notion of "conductor" , which, asymptotically, is the logarithm of the "usual" 
conductor used here.) 

Conjecture C. Suppose J 7 is a family of L-functions and, for F € T ' , let c(F) denote the 
conductor of F. With B = 1/2 for unitary families and B = 1 for symplectic and orthogonal 
families, we have 

max \F{\)\ = exp ((1 + o(l)) log £> log log £>) . (1.5) 

c(F)<D 

The implied constant depends only on T . 

For example, for the symplectic family of real primitive Dirichlet L-functions, L(s,Xd), where 
Xd — (.), we conjecture that 

max \L(\, X d)\ = exp ((1 + o(l)) ^/ log D \oglog d) . (1.6) 

Similarly, for the orthogonal family of Dirichlet series associated with holomorphic cusp forms, 
L(s,f), where / <E Sk(To(N)), the conductor is kN, so we conjecture that 

max |L(|,/)| = exp ((1 + o(l))V log L> log log L>) . (1.7) 

kN<D 

Note that Conjecture l£l contains Conjecture because any primitive L-function, L(s), has 
associated with it the unitary family 

T L :={F y (s) :=L(s + iy) | y e M} (1.8) 

with conductor c(F y ) ~ \y\. 

Our conjectures suggest that on the critical line the answer to the "1 or 2" question is "1" . 
Work of Montgomery and Vaughan ^7] and Granville and Soundararajan has suggested 

that the answer also is "1" on the 1-line. Thus in both cases, the maximum value the L-function 
attains appears to be closer to the fi-result than to the O-result. 

In Section [3 we use a rigorous approximation to the zeta-function due to Gonek, Hughes, 
and Keating [Hj to justify Conjecture El This approximation represents £(s) as a product 
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over primes times a product over zeros. We use characteristic polynomials of random unitary 
matrices to model the product over zeros, and a separate probabilistic model due to Granville 
and Soundararajan for the product over primes. The approximation to the zeta-function has a 
parameter which controls the relative contribution of the primes and the zeros. We show that the 
predicted maximal order of the zeta-function is the same independent of the choice of parameter. 
That is, whether we use only the primes, or only the zeros, or some combination of the two, we 
obtain Conjecture 1X1 

In Section we use recent conjectures for moments of L-functions to give an alternative 
justification for Conjecture^ Our approach also provides new limits on the range of validity of 
those conjectured moments. 

In Section 0] we modify the treatment in Section [21 to obtain Conjecture iBl 

In Section|31we describe how to extend our approach to obtain Conjecture^ We also describe 
some other approaches to obtaining these conjectures and then indicate possible arguments 
against the conjectures. 

Finally, in Appendix [X] we prove a theorem about random matrix polynomials that is used in 
Section 

We thank Andrew Granville and Soundararajan for allowing us to incorporate their work on 
extreme values using the primes, and we thank Andrew Booker, Brian Conrey, Hugh Montgomery 
and Doug Ulmer for helpful conversations. 



2. A PROBABILISTIC MODEL FOR THE ZETA-FUNCTION 

Gonek, Hughes, and Keating [H] have proved that if s = a + it, with < a < 1 and \t\ > 2, 
then for X > 2 and K any positive integer, 

C(s) = Px{ S )Z x {s) (l + O + °( X ~° lo S*)) . (2- 1 ) 

where 

*M'):=«p(E-^) , (2-2) 
\n<X 6 J 

A(n) is the von-Mangoldt function, and 

Zx(.s):=cxp^-^C/(( S -p)logA)^ . (2.3) 

Here the p are non-trivial zeros of £(s) and U(z) — u{x)E\ {z log x) dx, where E\(z) — 
I°° ^~w~ ^ w ^ s ^ ne exponential integral and u is any smooth function supported in [e 1-1 / x , e]. 

The parameter X controls the relative influence of the primes and the zeros. If X is large, there 
are many primes in Px(s), and only the zeros very close to s effect the product in Zx{s), while if 
X is small, the zeros further away from s make a contribution to Zx(s), but the number of primes 
in Px(s) is diminished. When X is not too large, we expect Zx and Fx to behave somewhat 
independently, and Gonek, Hughes, and Keating [H] give evidence of this. In the remainder of 
this section we describe probabilistic models for Px and Zx which, assuming independence, will 
give Conjecture [XJ In Section |2~TI we describe our model for the large values of \Zx\, establish 
some new results on the size of characteristic polynomials of random unitary matrices, and justify 
Conjecture ^ by choosing X small. In Section l2~2l we describe Granville and Soundararajan's 
model for Px and justify Conjecture El by choosing X large. Then, in Section l2~3l we combine 
Zx and Px, showing that intermediate values of X also lead to Conjecture lAl 
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2.1. A random matrix model for large values of Zx(h + it). Here we study the charac- 
teristic polynomial 

N 

Au(9) = det (/ - Ue- lB ) = JJ (l e 1 ^- 6 )) (2.4) 

n=l 

of a random unitary matrix U £ t/ (iV) chosen uniformly with respect to Haar measure. The 
characteristic polynomial Aj/(0) was first developed as a model for the Riemann zeta-function 
by Keating and Snaith m 0j building on it is argued that for t&T, Zx{\ +it), given 
by can be modeled by Au(0), where f7 € U(N) with 



AT = 



logT 



(2.5) 



eT logX_ 

We will prove a result about the value distribution and maximal size of |Aj/(0)| in Appendix lAl 
and use it to conjecture the distribution of large values of \Zx\- 

The largest value of |Aj/(#)| is 2 N , and values near this occur when U is in a small neigh- 
borhood of scalar multiples of the identity matrix. If X = e°( loglogT ) , this violates the known 
bound on |£(| + it)|, so our model for the large values of \Zx(\ + mus t do something more 
subtle than just take the maximum of |A[/| over all U £ U(N). 

If T and X arc thought of as fixed, then matrices of size N = f }°, sT v ] should model the 

O 7 L e 1 log X 1 

zeta-function as long as TjX e '' < t < T. If X > 2, say, then up to constants there are TlogT 
zeros in this interval. Therefore, in order to have the same number of eigenvalues, one needs 

M = [T log T/N] « exp (e 7 N log X) log X (2.6) 

matrices. Thus, one plausible guess for the maximum value of \Zx{h + for < t < T is 
K = K{M, N), where A" and M are given in l|2.5|l and l|2.6|l . respectively, and K is the smallest 
possible function of M and A^ such that 

^max^ max \A Uj (0)| < k\ -> 1 (2.7) 

as A^ — > co. Such a AT is found in Theorem l2.ll 

We have glossed over some issues here, but we argue that they are not essential. First, we 
have claimed that matrices of size N model Z x {\ + it) for T / X e ~' < t < T, whereas we want 
< t < T. However, if X ->■ oo, then [T/X e , T] will cover almost all of [0, T], and so should 
capture the maximum. Secondly, we have been slightly cavalier in dropping the condition that 
A^ should be an integer, which will have an effect on the number of matrices, M, we maximize 
over. However, as we will see below, the answer depends only on the logarithmic size of M, so 
this is not a serious problem. Finally, we remark that the placement of e 7 in our definitions of 
A" and M is actually irrelevant: our heuristics are sufficiently robust that increasing A^ by any 
fixed constant and decreasing M correspondingly leads to the same conjectured maximum. We 
include the e 7 factor to be consistent with jS], where the precise choice of A" does matter. 

We now find an explicit K satisfying (|2.7|l . 
Theorem 2.1. Fix S > 0. Let M = exp(A^), with 6 < (3 < 2 - 5, and set 



K e {N) = exp 



( (^1 - \p + V log A/ V log Af) . (2.8) 



If U\, . . . , Um ar £ chosen independently from U(N), then as N — > oo, 



max max |A[/. (#)| < KJN) > — * 1 (2.9) 
1<j<M e 3 1 



for all e > and for no e < . 
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Proof. Note that by the independence of the Uj, 

^ax^ max ^(61)1 < if e (JV)| = P {max \A V (9)\ < K e (N)^ (2.10) 

and, for this to tend to 1 as TV — ► oo, we must have 

MlogTl-P jmax|A,7(0)| > K e (N)^ -> . (2.11) 

Thus, the proof of the theorem (and all similar ones in this paper) requires knowledge of the 
tails of the distribution, and this is given by Lemma fA. II If M — exp(N ,s ) and K — K £ (N) is 
as in (|2.8|l with < (3 < 2 and e > — yl — (3/2, then by Lemma [A. II one easily finds that 

M log (l - P {max ^(9)] > K £ (N))) = exp(T^) exp + o( l)) 

(2.12) 

-> 

as TV — ► oo for all e > 0, but for no e < 0. □ 

To summarize, we use the characteristic polynomials Ajj(9) of random unitary matrices U G 
U(N) to model Z x {\ + it)- To model the large values of \Z X {\ + \t)\ for t € [0, T] we choose TV 
as in (|2.5|l with logX < (log N) A for some A, and we take about 

M = N c exp (e 7 /V log X) (2.13) 

different matrices. Here c > is fixed, and we include it to allay concerns that choosing too 
few matrices may miss some large values. With these values for M and TV, it follows that (3 in 
Theorem 12 .11 is ~ 1, and this leads to the following conjecture: 

Conjecture D. If 2 < X < log A T, then 



max \Z X {\ + it) | = exp (1 + o(l)) \ \ logTloglogT (2.14) 
te[o,T] z \ v £ J 

as T — > oo. 

We can now complete our argument for Conjecture 1X1 
Justification of Conjectured^ By the prime number theorem and (|2.2|l . we see that 



Thus, if X = O(logT) and T/X fiT <t<T, then 

M1 + -)-o(«p(o^)). (2-16) 

Combining this, (|2.1|1 . and Conjecture ^] we obtain Conjecture ^ □ 

The argument above essentially splits the critical line into blocks of size 1, maximizes over 
each block, and then finds the maximum of the maxima. However, one might instead wish to 
sample the critical line at many evenly spaced points. If they are not too sparse, then a value 
close to the global maximum in [0, T] will be found. The following lemma explains why this is 
the case. 

Lemma 2.2. Suppose +i*o)| = m T '■= m ' AX te[a.T] \C(h There is an absolute constant 

A>0 such that if \t - t \ < A/logT, then + it)\ > ~m T . 
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Proof. We can estimate the size of the derivative of the zeta- function near | + iio using Cauchy's 
theorem. If we integrate around a circle of size 1/logT and use the functional equation, we 
find that there exists an absolute constant ci such that if \s — (4 + iio)| < Ci/logT, then 
C'(s) <gC wit log T. This gives the lemma. □ 

As a random matrix model for \Zx{\ + ii ) | when it is sampled at evenly spaced points, one 
might consider the largest value of K = K(N, X) such that 

p{ max \A V (0)1 < k\ 1. (2.17) 

[l<j<V<= exp(AflogX) 3 J 

The following theorem determines K explicitly as a function of N and X and shows that such 
sampling is sufficient to capture the large values. 



Theorem 2.3. Let M = N°e logA , where 2 < X < N and c> is fixed. If 



K E (N) = exp ^(-L + e)^/NlogNlogX^ 



(2.18) 



and U±, . . . , Um are chosen independently from U(N), then as N — > oo, 



^JAuM <K e (N)\->l (2.19) 



/or all e > and /or no £ < . 

Proof. As in the proof of Theorem 12. II since the Uj are independent we have 



log I 



max |A^(0)| < K E (N)\ = Mlog(l -P{|At,(0)| > K e (N)}) . (2.20) 

1<J<M J 

Theorem 3.5 of H2 asserts that if 5 > is fixed and exp(A 5 ) < K < exp(A 1 -' 5 ), then 

F{\A um > K} = exp (- ^-^— (1 + 0,(1))) (2.21) 

as N -> oo. Hence, if M = A c exp(AlogA) and K e (N) is given by J^TTSJ, then the left-hand 
side of H2.20f> tends to zero for all e > 0, but for no e < 0. □ 

Note that the statements of Theorems 12.11 and 12 . 31 are almost identical and, in particular, one 
can capture the largest values of |A[/|, hence \Zx\, just by sampling at individual points; it is not 
necessary to find the maxima of the individual polynomials. This is significant for our modeling 
of the prime contribution Px , for in that case we are only able to sample at individual points, 
and there is nothing comparable to a sequence of polynomials over which we can maximize 
individually. 

2.2. Probabilistic model for large values of Px- The material in this section was provided 
to us by Granville and Soundararajan. 



First note that 



= E ^ + 0(loglogA). (2.22) 



p<X v 



Hence 



Px{\ + it) = exp £ -4^ x exp(0(log log A)) 
\p<xP 2 J 
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= exp(P£(i + rf)) x ex P (0(loglogX)), (2.23) 
say. We will see that this approximation is adequate as long as X = exp(o(-y/ log T log log T) ) . 

The method is based on treating the p~ lf as independent random variables. The large values of 
Pxih+rt) can then be obtained from the following lemma, the proof of which involves calculating 
the moments of the distribution. 

Lemma 2.4. Let {zj} be a sequence of independent random variables distributed uniformly on 
the unit circle and let {aj} be a sequence of bounded real numbers such that for all n > 3, 

-L £ ay-0 (2.24) 

as J — > co, where 



V J l<j<J 



Vj:= 4 ( 2 - 25 ) 

i<j<J 



Then, as J — > co, the distribution of 



Yj := me a i z i ( 2 - 26 ) 

i<j<J 

tends to a Gaussian with mean and variance \Vj. 

Applying the lemma to Px{\ + rt) with Zj = pj 1 * , where pj is the jth prime, and aj — 1/ \fp~j, 

we see that Vj ~ loglogX. For X = exp(^/ log T) we model the maximum of \P X {\ + d)\ by 
independently choosing Tlog c T values of t. By (|2.23H this yields 

4 max \P X {\ + it)\ = exp ^(1 + o{\))^\ log T log log Tj . (2.27) 



If X = exp(ydogT), the method of Section I2~T1 predicts that 



max \Z x (±+it)\ = 0(cxp( v /logr)) . (2.28) 
These estimates together with ()2. ID give an alternative justification of Conjecture lAl 

2.3. Combining Zx and Px- If X = cxp(log" T) with < a < |, then the largest values of 
\Zx | and \Px\ are approximately the same size and both will contribute to the largest values of 
|£(i +i£)|. Specifically, applying Theorem 12 . II with N = logT/logX (so that Z x is modeled by 
A[j) and M — Tlog X (so that we sample enough characteristic polynomials to cover the critical 
line between t = and t = T), the previous analysis using characteristic polynomials predicts 
that \Zx{\ + W)| gets as large as 

exp ^-^(1-20) log Tlog log 7^ , (2.29) 
and \Px(h + S e ^ s as large as 



exp (ya log T log log Tj . (2.30) 

The product of these two quantities is larger than our conjectured maximum of |C(| + 
as it should be, because we do not expect \Zx\ and \Px\ to attain their maximum values 
simultaneously. Instead of multiplying the maxima, we must find the distribution of the large 
values of the product \ZxPx\ in order to check that our method is consistent throughout the 
range < a < \. This is a calculation involving the tails of the distributions of Zx and Px ■ 
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In Section \'2.2\ we saw that, conjecturally, if X = exp(log Q T), then the tail of logPx has 
the same distribution as a Gaussian with mean and variance <jp — \a log log T. For Zx, 
Lemma [A. II states that if 8 > is fixed and 8 < \ < I — 5, then 

max|A a (g)| >exp(^ A )| =exp(^- (i ^ Qgjv (l + (l))j . (2.31) 

Since Ajj models Zx, the lemma provides the tail of the distribution of log j^xd+ii)]. Assuming 
that Zx and Px are essentially statistically independent, we should expect the distribution of 
log \ Zx \ +log \ Px\ to be the convolution of the two distributions. Hence, for the tail we convolve 
the tails of the two distributions. Thus, for large K we expect that 

imeas{0<t<T : log |P x (i + if)| + log \Z X (3 + it)| > logif} 

CXP {- { a log log T + (l- Q )loglogT-log, J (1 + J ^ 



— OO 



e X p(-f K (x)(l + o(l)))dx, (2.32) 

say. By the saddle point method, this equals 

exp(-f K (x )(l + o(l))), (2.33) 
where xq is such that f' K (xo) — 0. If we write K = exp (d^/ log T log log T^j , then solving 

2 (x - d^logTloglogT^j 2xq 



= f' K (xo) = 



a log log T (1 — a) log log T — \ogxo 

x 



((l-^loglogT-log^o) 2 



(2.34) 



yields the solution xq ~ d(l — 2a) a/ log T log log T (which is justified so long as < a < 1/2). 
Thus, 

f K (x ) = (2 + o(l)) d 2 logT, (2.35) 
and this leads to the following conjecture: 

Conjecture E. For d > fixed and T — > 00, we have 

imeasjo < i < T : |C(j + i<)| > exp^^logTloglogT)} = exp (-2d 2 logT(l + o{\))) . 

(2.36) 

By Lemma l2~2l \((^ + it)\ is close to its maximum value over a window of size C/logT, so 
we wish to find the smallest d such that 

meas{0<i<T : + it)\ > cxp^v/logTloglogT)} < — L;, (2.37) 



that is, such that TlogTexp(-2d 2 (l + o(l))logT) < 1. This happens if d = y ^ + £ for any 
e > 0, but for no e < 0. Once more this gives Conjecture El an d this time in a way that is 
independent of the choice of X = exp(log a T) for < a < 1/2. 

3. Bounds based on conjectures for moments 

In this section we obtain ConjectureElby using conjectures for moments of the zeta-function. 
Our method also leads to limits on the possible uniformity of the conjectured moments. 

Our approach here is based on the work of Conrey and Gonek |H] . Let 

rriT ■= max (3.1) 
te[o,T] 
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and note that we have the trivial inequality 

m| fc > ^ T |C(i + it)| 2fc d^. (3.2) 

It follows that estimates for the right-hand side of the inequality imply lower bounds for the 
maximum size of the zeta-function. 

Keating and Snaith |13| used random matrix theory to conjecture that if k > —1/2 is fixed, 
then as T — ► oo, 

' ^ 1C(| + i*)l 2fc dt ~ ^ + ]\ a(k) log fe2 T, (3.3) 



T J Q ISV2 n G(2fc + 1 

where G is the Barnes G-function, and 



p 

prime 



This conjecture is for fe fixed, but we would like to let k — > oo, because the 2fcth root of the 
right-hand side of (|3.2|l then actually tends to my. Thus, we would like to know how large k 
can be as a function of T. 



Conrey and Gonek showed that if formula 1)3. 3|) holds for k = log T / log log T then 



log log T 

and if it holds for k as large as log Tj log log T then 



mr>exp(C 1 J 1 l0 f T ^ ) , (3.5) 



m T > exp ( C 2 l0gr ) , (3.6) 



' log log T 

where C\ and C2 are given explicitly. Hughes gave a convexity argument to show that 
formula (|3.3II must fail before k = log T/ log log T. However, using the last point at which 
convexity holds for l|3.3|l . one still obtains l|3.6|l . but with a smaller constant Ci- 

In all of these cases one only requires a lower bound for the right-hand side of l|3.3|l . Our 
approach here is to use the mean value formula l|3.3|l to obtain upper bounds instead of lower 
bounds for my. As a consequence, we also obtain restrictions on the possible range of validity 
of (|3.3|l for k growing with T. Specifically we prove the following: 



Theorem 3.1. Formula (I3.3|l does not hold for 



(3.7) 



/or any fixed e > 0. 



Our method allows us to get upper bounds for tut and, in particular, we obtain 
Theorem 3.2. If formula holds for k = log 5 T for some S < i, then 

m T < exp (log 1 " 5 r) . (3.8) 



Moreover, if formula (|3.3|) holds for k = yj 2 log T / log log T, then 



mT « exp ^T W+ ^y , (3.9, 
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and if formula (I3.3|l holds for k = yj 8 log T / log log T, then 



»>, -M- ( y 4 logTloglogT + Q (VWWW**T\\ 

\ y/\0g\0gT II 



Theorem 13.21 says that if formula (|3.3|l holds true until k — */ 8 log T/ log log T (after which 
we know it must fail), then we have 

log log log T 



m T = exp (y I logTloglogT \l + ( J j j (3.11) 

Note that this implies the conjecture found in the previous section. This is not surprising 
because, as we will show, the arithmetic factor in the conjectured moments is smaller than the 
other factors. Thus, this bound is coming just from the random matrix model. 

If the true order of the zeta-function is larger than the bound in lift , then one would like 
to know where our calculation fails. Since formula l|3.3|l is only the leading order term in the 
asymptotic expansion for the 2Mh moment of the zeta-function, it is possible that the lower 
order terms dominate when k — log 5 T. However, this is unlikely. In the random matrix case, 
the 2fcth moment is given by 

F/lA (ft\ 2k \ G 2 (k + l)G(N + l)G(N + 2k + l) 

M 1 ^ 1 J = G(2fc+l) GHN + k+1) (3 - 12) 

_G 2 (fc + l) / _ Ar , fc 3 7k 



G (2k + l)^{ k ^ N+ N-T2N^ + 2m + ---)' (3 ' 13) 

and one sees that the first term dominates even for k as large as N S T provided that 8 < 1. 

In the zeta-function case, the complete main term of the 2fcth moment has been conjectured 
(see [S], Conjecture 1.5.1). One can check that the contribution from the primes is bounded 
by exp(cfc 2 ), which is insufficient to affect the estimate for rar- Thus, if our conjecture for the 
growth of + \t)\ is incorrect, then the main term in the mean value must take a new form 
for k = log 5 T for some 8 < h. If Ijl.lfl is the true maximal size, then by equation (|3.2U|) the 
conjectured mean value can only hold for k -C log log T. 

Our main tool for finding upper bounds is the following lemma. 
Lemma 3.3. For all positive real k we have 

mf <2 2fc logT [ |C(i+it)| 2fc dt, (3.14) 
Jo 

where the implied constant is absolute. 

Proof. Suppose that + ito)| = %i where < to < T. By Lemma \'2. 21 there is an absolute 
constant A > such that if \t - t \ < A/\ogT then |£(| + it)\ > \m T . This gives 

-io+A/logT 

a 2 

to-A/logT 

>- 

as claimed. □ 

Proof of Theorems \3 . 1\ and \3. 6 A By Lemma [3.31 and (|3.2|l there exists an absolute constant C 
such that 

l±J o m+it)\ 2k dt\ <m T <2(CTlogT) 1 /2Ml^ | C (i +lt) |^ d ij . (3 . 16) 
We use these inequalities to prove Theorem 13 . 21 first . 



\Ca + it)\ 2k dt > / |c(| + rt)| 2fc d* 
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The asymptotic expansion for the Barnes G- function (see Barnes fSj) gives 
(/( ' ' ' - exp (k 2 (-\ogk+l-2 log 2^ - -L log + 1 log 2 + C'(-l) + O ( ' 



(3.17) 



G(l + 2Jfe) V X \ b 2 " y 12 " 12 " * v ' V*, 

for fc > 1. Furthermore, Conrey and Gonek [2] have shown that 

loga(fc) fc 2 log(2e 7 logfc) + o(fc 2 ) for fc -> 00. (3.18) 

Thus, if holds, then 



( t \ 1/21 

log f ^ |C(§+i*)f di) =i£loglo g T-^log£ + 0(£loglog£). 



(3.19) 



It therefore follows from l|3.16[l that 



log T + log log T 1 1 
logm r < S ^ B S +-£loglogT--llogl + 0(lloglogl). (3.20) 



Setting ^ = log 5 T, we find that 



//// ( ixi» ( ilog^T + OQog^TloglogT) j . (8.2:1 ) 



which gives the first inequality in Theorem 13. 21 Setting I — cJ ^°fj ' T , we find that 



which is minimized by taking c = y^2. This gives the second estimate in Theorem 13.21 
By (|3~T^|) and (|3~T3|) we also have 

m T > exp QfcloglogT- ifclogfc + O(fcloglogfc) j . (3.23) 



If we set k = c\l , lo f T „ , we find that 

1 ' log log T 5 



v/logTlogloglogT 



mT » exp ^VlogTloglogT + O ( J 7==^ J J • ( 3 - 24 ) 

Choosing c = 2-^/2, we obtain the third inequality in Theorem 13.21 If c > 2-^/2, this contra- 
dicts ()3.22(l and thereby establishes Theorem 13. II □ 



4. Bounds for S(t) 

Recall that S(t) is the error term in the counting function for the number of non-trivial 
zeros of the zeta-function with imaginary part less than t. It may also be expressed as S(t) = 
iamlogC(^ + it) (see Titchmarsh |T§J). Since ({^+it) is essentially P x (^+it)Z x and Z x 
is modeled by Ajj(0), one would expect that if X is sufficiently small, so that the contribution 
from Px is negligible, then i?mlogC(^ + it) too can be modeled by random matrix theory, in 
particular, by 

-3m\ogAu{0). (4.1) 
Evidence for this is presented in ^3]. This is the basis for our Conjecture iBl 
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Theorem 4.1. Set 

(4.2) 



K e (N) = ^-L + £ ^J y/NlogN 



and M = N c e N , where c > is fixed. If Ui, . . . , Um ar ^ chosen independently from U(N), then 
as N — > oo, 

max 3m log Ay . (0) < KJN) \ -> 1 (4.3) 

l<j<M J 

for all s > and /or no £ < . 

Proof. This follows along the same lines as the proof of Theorem 12.31 Since we are making 
independent choices, 

max Jm log Ay . (0) < K \ = P {3m log Ay (0) < K} M . (4.4) 
1<3<M J 

For this to tend to 1, we need 

MlogP{amlogA c/ (Q)<^}-*-0. (4.5) 

By Theorem 3.6 of Hughes, Keating and O'Connell if if = N x , where 6 < A < 1 — S and 
S > is fixed, then 

PpmlogA^O) <K)} = 1- exp (- ^ R (1 + 0ff (l))) . (4.6) 
One can easily check that if K £ (N) is given by l|4.2() . then l|4.5() holds for all e > 0, but for no 

£ < 0. □ 

Conjecture^] now follows in the same manner as the justification of Conjecture^ that is, by 
controlling the prime contribution from 2m\ogPx{\ + it). 



5. Other families and other arguments 

5.1. Other families: symplectic and orthogonal. The analogue of Gonek, Hughes, and 
Keating's approximation to the zeta-function has not yet been extended to the case of other 
L- functions near the critical point. However, it is believed that the characteristic polynomials 
of symplectic (or orthogonal) matrices model the central value of L-functions taken from a 
symplectic (or orthogonal) family of L- functions ^31 • Thus, the methods developed in section l2~T1 
can be applied. Moreover, we can still estimate the maximal size of critical values by using a 
partial Euler product and modifying the method of Granville and Soundararajan. Finally, we 
can also apply the method involving mean values. 

We give as an example finding the large values of the characteristic polynomials of the sym- 
plectic group at the critical point. The orthogonal family is treated in an almost identical way. 
The characteristic polynomial of an N x N symplectic matrix (iV must be even) with eigenvalues 

e ±i6 n ig 

N/2 

Z(U,0) = l[(l-e^)(l-e-^). (5.1) 

3=1 

Keating and Snaith XI] calculated the moment generating function and found that 

f <7(Tj frt'i - ? N ° TT T{N/2 ± 3 + 1)r( ' s ± I ± 1/2) an 

E Sp{N) {Z(U,0)}-2 lJL r(i + 1/2)r(s + iv/2+i + ir (5-2) 
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A long but straightforward calculation using Stirling's asymptotic series for the gamma function 

(l-A)logAT' 



shows that if S > is fixed and 5 < A < 1 - S, then for A(N) = N x , B(N) = ( .l N , and 



fixed s > 0, we have 

lim l\ O gE Sp{N) {z(U,0y B / A } = l S 2 . ( 5 - 3 ) 
TV — >oo / j \. J z 

From this, large deviation theory (for example, see [7]) allows us to deduce that if exp(N s ) < 
K < expiN 1 - 5 ), then 

F Sp{N) { Z(U, 0)>K} = exp (- 2l0gJV l0 _ g ^ gl0gj , (l + o,(l))) (5.4) 

as N — ► oo. Comparing this with (|2.21(l . the analogous statement for the unitary group, we 
note the extra factor of 2 in the denominator. This difference explains why the constant B in 
Conjecture ICl equals 1 rather than 1/2. 

We now see, by methods identical to those of the previous section, that if M = N c e N for any 
fixed c > 0, and if K e = exp nl + s)\/N log N^j , then if Ui, . . . , Um are chosen independently 
from Sp(N), 

P Sp{N) < K e (N)} -» 1 (5.5) 

as N — > oo for all e > and for no e < 0. 

Consider for instance the family of all quadratic Dirichlet L-functions L(s, Xd)- For characters 
with modulus around D, random matrix theory suggests (see Keating and Snaith ^31) that 
N = log D is the correct identification between the size of the matrix and the conductor (though 
note that in ^3] N is half the size of the symplectic matrix). Furthermore, it is well known that 
there are asymptotically 6D/ir 2 primitive discriminants less than D. Thus, we conjecture that 



max =exp((l + (l))Vlog£>loglogD). (5.6) 

Xd real 

Similarly, the moment generating function has been calculated for the orthogonal case (see 
^31) an d ; if N is even, we have 

E SO[N) {Z(U,0)} = 2 [[ T{j _ 1/ms + N/2 + j _ iy (5.7) 

Equations (j5.3|l and l|5.4|l apply to the orthogonal case without change, so by the same reasoning 
as previously, if M = N c e N for any fixed c > 0, and K e — exp f(l + e)y/N log N^j , then 

P so(N) ( max .Z(U jt 0) < K £ (N)\ -» f (5.8) 

as N — ► oo for all e > 0, but for no e < 0. 

Next we consider how to adapt the Granville-Soundararajan argument involving the product 
over primes to the symplectic case. We require the following lemma. 

Lemma 5.1. Let {xj} be a sequence of independent real random variables with mean and 
variance 1, and let {a., } be a bounded sequence of real numbers such that for all n > 3, 



as J — > oo, where 



Vj x<j<J 



v.,:= a r ( 5 - 10 ) 
l<j<J 
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Then as J — > oo, the distribution of 

Yj:= a J x i ( 5 - n ) 

i<j<j 

tends to a Gaussian with mean and variance Vj. 

Just as in the treatment involving characteristic polynomials, the fact that the variance is 
Vj for these families instead of Vj/2 leads to the constant B — 1 in Conjecture [UJ instead of 
B = 1/2 for the unitary family dealt with previously. 

5.2. Other arguments, for. Although the conjectures in this paper are based on very recent 
work, Hugh Montgomery has pointed out to us that a similar conjecture can be obtained by 
viewing log |C(5 + 1 ^)I as a Gaussian distributed random variable with variance CloglogT, where 
one estimates my by sampling T A times. 

Soundararajan suggests a different way to use the moments of the zeta-function to conjecture 
an upper bound. The proof of Lemma 13.31 showed that if 

imcas{<e [0,T] : |C(| + if)| > r} < ± (5.12) 

then rriT < 2t. For when + it) | is very large, it must remain large over an interval of size 
c/logT. Now 

T 2k me&s{t e [0,T] : |£(± + it)\ > r} < [ T + it)\ 2k dt , (5.13) 

Jo 

so if equation (|3.3|l holds, then we have 

imcas{tG[0,T] : |C(| + it)\ > r} < r~ 2fc ± ^ (log Tf . (5.14) 

The right-hand side is less than c/{T\ogT) (which means there is only one place where the 
maximum occurs) when 

r^expf^ + ^loglogT-^logfcV (5.15) 



The minimum of this is exp [\J j log T log log T^j and it occurs when k = \J2 log T j log log T. 

5.3. Other arguments, against. We now discuss potential arguments against the conjectures 
in this paper. One possibility, so fundamental that it cannot be addressed, is that the large 
values of an L-function may be so rare that these statistical models cannot detect them. Indeed, 
since these are problems in number theory, there may be number-theoretic constructions of large 
values which contradict our conjectures. The two examples below, due to Brian Conrey, suggest 
the kinds of things we have in mind. 

The first argument invokes an analogy with the divisor function d(n) = Y^,M n 1 and the related 
function cu(n) — 2p| n 1' where here the sum is over prime divisors of n. Since uj(n) is log log n 
on average and has a Gaussian distribution, the relation d(n) — for square-free n might 

lead one to conjecture that for such n, d(n) is bounded by 



exp (c \J log nj log log n) . (5.16) 

However, we know how to construct large values of d(n) , and it is easy to see that for n squarefree, 
d(n) can get as large as 

exp(clogn/ log log n). (5-17) 

The second argument concerns the Fourier coefficients a n of cusp forms. For integer weight 
cusp forms, rescaled so that for p prime we have \a p \ < 2, the coefficients a n can get as large as 

exp(clogn/ log log n). (5.18) 
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In other words, they can get about as large as d(n). The question is: can the coefficients of half- 
integral weight forms also get this large? If they can, then our conjecture on the maximal size of 
the critical values of a symplectic family of L-functions is incorrect. For if / € Sk(To(N)), then 
Lf(h,Xd) = c d/\/^' w here Cd is a Fourier coefficient of the half- integral weight form associated 
with / by the Shimura correspondence. 

Our methods cannot be directly applied to produce a version of Conjecture El for symplectic 
and orthogonal families. But if one were to believe that for a family T of L-functions with c(F) 
denoting the conductor of F € 

3mlogF(i) 



limsup K1 = y/B (5.19) 

c(F)^oo V log c(F) log log c(F) 

where B = 1/2 for unitary families and B — 1 for symplectic and orthogonal families, then since 
the rank of an elliptic curve is related to the order of vanishing of its associated L-function, this 
could lead to new information about large ranks. That is, if (|5.19|) is true, then it suggests that 
for rational elliptic curves we have 

rank(E') 

limsup V ' = 1, (5.20) 

c E ^oo y/ log C E log log C E 

where ce is the conductor of E. Note that this is smaller than the ranks of elliptic curves found 
by Ulmer [2U| in the function field case. 

Appendix A. The tail of the distribution of maxg |A[/(0)| 

Here we prove the random matrix polynomial result used in Section 
Lemma A.l. If 5 > is fixed and 5 < A < 1 — 8, then 

maxIA^)! > exp(A A )| = exp (- —J^L—^ Q. + (1))\ . (A.l) 

Proof of Lemma \A.l\ Bernstein's inequality for polynomials implies that for any matrix U, 

max | A'^ < A max | Ay (6>) | . (A.2) 

9 9 

Thus, if is a point at which the maximum of |A[/(0)| occurs, and if we indicate the maximum 
by mu, then for \6 - <j>\ < 1/N, 

\Au{9)\ >m v - \8-<b\N mu . (A.3) 



It follows that 



,1/AT 

\^u{d)\ 2k d9 > / (1 - \x\N) 2k dx (A.4) 

J-l/N 

= m $] L ^_ (A5) 

u N 2k + 1 y ' 
Combining this with the trivial lower bound for mu, we find that 

— / \Au(6)\ 2k d6 < < ~^-N / \A v (9)\ 2k dO. (A.6) 
2tt Jo 2 J 

This bound holds for any matrix. We now average over all N x N unitary matrices with respect 

to Haar measure. That is, we calculate the expectation E^r of \A(j(0)\ 2k - Set 

E N {\Au(6)\ 2k } = M N (2k). (A.7) 

Keating and Snaith |13| have shown that 

M (2k) G 2 (k + l)G(l + N)G(l+N + 2k) 

MNi2k) - G (2k + l) G 2 (l + N + k) ' (A - 8) 



16 



DAVID W. FARMER, S.M. GONEK, AND CP. HUGHES 



where G is the Barnes G-function. Note that this is independent of 9. Therefore, by (|A.6J) 

M N (2k) <E{ml k } < n{2k + l)NM N (2k). (A.9) 

Hughes, Keating and O'Connell [Hj have shown that if A(N) = N x with 6 < A < 1 — S and 
S > fixed, and if 

B ^ = (l-A)logiv ' (A.10) 



then for s > 0, 



1 ( sB i N )\ _ 1 „2 



^BiN) l0gMN {M^))^r- ^ 



Since 

*logM N (sB(N)/A(N)) < ^—^ogE{mf (N)/AiN) } 



~ B^b6M N (8B(N)/A(N)) +0 (^|^ ) ^ (A. 12. 



we conclude that for s > 0, 



1 , W J ^i3(A)log(max 9 |A^(0)|) ^ 1 2 

v 1 ™ S(A0 l08E ( eXp { AjN) J} = I S • (A ' 13) 

From this, large deviation theory (see, for example, \7\) allows us to deduce that 

lim — — logp{logmax|Ay(0)| > A(N) \ = -1. (A.14) 



N^oo B(N) 

Inserting A(N) = N x and B(N) fr om 1A.10I). we obtain the statement in the lemma. □ 
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